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1 Introduction 

In this paper we consider a class of scalar IVP's for linear dispersive equa- 
tions with fast temporal and spatial scales subject to highly oscillating 
initial data. The Cauchy problem of the Schrodinger equation serves as 
a typical example 

iedt^^ = -—Aip' + V{x)i;'' = 0, xeR'^,teR (1.1) 

^"(x,0) = Ai{x)e'^'^'-''>/', xeR'^ (1.2) 

where e ~ (the scaled Planck's constant). The small parameter e rep- 
resents the fast space and time scales introduced in (1.1), as well as the 
typical wave length of oscillations of the initial data. We are interested in 
the high frequency limit of these equations, which is usually referred to 
as "geometrical-optics" . In the special case of the Schrodinger equation 
with vanishing Planck's constant this is precisely the " (semi-) classical 
limit". It is well known that the considered equations propagate oscil- 
lations of wave lengths e which inhibit from converging strongly in a 
suitable sense. 

Thus the short- wavelength-asymptotics e — > is by no means straightfor- 
ward, in particular since the physical quantities of interest (observables) 
are quadratic in ip^ . 

The usual way to tackle the problem is the geometrical optics - or WKB- 
Ansatz (Wentzel-Kramers-Brillouin, see [[Ke|), which consists of repre- 
senting the solution in the form 

ip^x, t) ~ A^x, t) exp (^^S{x, t)^ (1.3) 

where A^ and S are realvalued, > and in general A^ ^ A -\- eA\ -|- 

£2^2 + .... 

Then after inserting the above representation into the equation and by 
considering, as a first approximation, only the lowest order terms, one 
finds that: 

• the phase 5 is a solution of a nonlinear first order equation of 
Hamilton- J acohi type 

• the (zeroth order) amplitude satisfies a linear first order PDE (called 
transport equation) that can be brought into the form of a conser- 
vation law for the energy density n = A"^. 

A severe drawback of this method should be noted. The obtained non- 
linear equations do not have global, i.e. for all t G M, smooth solutions 
(except for some special initial data). In other words the system in gen- 
eral develops singularities in some finite time tc {"break time"). The 
formal expansion method clearly can only be justified for smooth, i.e. 
sufficiently often differentiable, functions S and and thus the Ansatz 
( |1.3D breaks down at points where the first singularities occur. These 
singularities are called focal points, or more generally caustics, since, as 
we will see, the energy of the wave becomes infinite there. 
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A alternative point of view on this problem is given by considering the 
so called Quantum Hydrodynamic System, which is obtained by plugging 
into equation (1.1) the ansatz (1.3), with S = S^. Defining := (A^)"^ 
and := n'^VS^, one gets (after seperating real and imaginary parts) 



0, 



dtf + div 



f ® f 



A 



:i.4) 
:i.5) 



This system is exact, i.e. equivalent to the Schrodinger equation (|l.l| ), 
c.f. | GaMa | and well posed for all t G M due to the third order disper- 
sive regularization term. However for e = 0, where the system simplifies 
to the zero temperature Euler equations, singularities appear and the 
equations cannot be used after them to describe the propagation of the 
energy density in the high frequency limit. (This approach is used in 
particular in one space variable for the non linear Schrodinger equation, 
i.e. = IV'^P, see e.g. [ |LaL(^ ].) 

A natural alternative to the standard WKB-method is seeking multival- 
ued phases corresponding to crossing waves. This means that in general 
for every fixed (x, t) G M"^ x R, which is not on the caustic, one tries to 
construct a (maybe infinite) set of phase functions {Si{x,t)}, i G / C N, 
each of which is a solution of the Hamilton-Jacobi equation in a neigh- 
borhood of (x,t). This set is referred to as the multivalued solution of 
the Hamilton- J acohi equation and it induces the multivalued solution 
{nj(x,t)} of the conservation law for the energy density. 
Historically this problem was studied by P. Lax, D. Ludwig, V. Maslov, 
J. Duistermaat and others (see |Dij], |Ki], |La], |Lu], | Masl| ]), who showed 
that Fourier- Integral operators and Lagrangian manifolds in phase space 
provide a uniform description of of the behavior of il^^ . (For applications 
in semi-classical quantum mechanics see the expository article of Robert 
[p!c|| .) The qualitative study of the multivalued solution is accomplished 
using geometri cal te ch niqu e s of si ngularity theory and contact geometry 
(see e.g. @, jAVGl] , pu[, f&ISl ). 

A considerable amount of work has been done in recent years on con - 
struc ti ng num erically the multivalued phase function (see, e.g., | Bel ], 
pe2i , pKl^ , gini). We will not cover the arising numerical 

questions in this paper, instead we refer the interested reader to these 
references. 

In the last decade the use of Wigner functions and Wigner measures 
has drawn increasing interest, in particular its application to the semi- 
classical hmit of Schrodinger equations ( pPa| , |MaMa|] , [[MMP| , |MPP|| ) 
and the homogenization of energy densities of dispersive equations (e.g. 
[|BCKP[| , [ paMa| , |Ge| , pMMP| ), mostly by groups in Europe. Indepen- 
dently, groups in the US used Wigner functions, too, wi th some e mphasi s 
on waves in random media and applied problems, e.g. [BKPR|, [PaRh]. 
The Wigner transformation provides a phase space description of the 
equations of the problem, which is extremely useful for the asymptotics 
since it "unfolds" the caustics (de-projection in phase space). Another 
advantage is that the high frequency limit, using Wigner functions, needs 
much lower regularity assumptions on Aj and Sj than the (generalized) 
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WKB-method. This is not only of purely academic interest since very 
often in concrete physical models Cq^ initial data are simply not avail- 
able. The price paid for the analytical convenience lies in the doubling 
of the dimension, i.e. the Wigner function is defined on M^"^. 
This work studies the connection between the WKB-method and the 
Wigner transformation or, in other words, represents an alternative ap- 
proach to WKB-asymptotics. It is organized as follows: 

• Section 2 is devoted to the setting of the problem. There we also 
give a short review of the traditional (naive) WKB-method and its 
generalization using FIO's. 

• In section 3 we present the main theorems on Wigner transforms 
and show how they can be used to obtain a semi-classical phase- 
space description. 

• This limiting phase-space regime is analyzed in section 4, which is 
the most important part of this paper. 

• Examples are studied in section 5 to illustrate the results of the 
foregoing section. 

• Finally in the Appendix in section 7 we give an example with non- 
global Hamiltonian flow and comment on the arising fluid type equa- 
tions, which generalize the zero temperature Euler equations. 

2 Setting of the problem and the WKB-method 
2.1 The model equation 

We consider the following initial value problem (generalized linear dis- 
persive model) for an anti-selfadjoint scalar pseudo-differential operator 

sdtip^ + iH^{x, sD)^p^ = 0, xeR'^,teR (2.1) 
subject to the highly oscillatory (WKB) initial data 

ijj'ix, 0) = ^/^^e'^'^'''>/', xeR'^ (2.2) 

where ip'^{t,x) is a scalar L^-function on M*^, D := — iV^, e G (0,eo] a-nd 
H{x,eD)^ is the Weyl-operator associated to its symbol H{x,e£,) by 
Weyl's quantization rule: 

Definition 2.1. 

H'^{x,sDMx):=-^ [ [ H(^,eA^{y)e'^^-y)^d^dy 

(2.3) 

The convenience of the Weyl-calculus lies in the fact that a scalar Weyl- 
operator is formally selfadjoint iff it has a real-valued symbol. 
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We have chosen this particular form of pseudo differential-calculus in or- 
der to be consistent with the usual framework of the Wigner-functions 
introduced in |GMMP| ]. Note that in case H is a sum of separate terms in 



X and the "Weyl symbol" and the "left symbol" of the classical Fourier 
multiplier coincide. 

Remarks. 

• The general framework of pseudo-differential operators allows us 
to include also non-local Hamiltonians, like the one appearing in 
example (iii) below, in our discussion. 



• Further note that the time and spatial scales of ( |2.1| ) are " fast" , since 
the small parameter e multiplies the time and spatial derivatives. 

We shall use in this text the following definition of the Fourier transform 
T : S(R'^) — > S( 



fiO := := / /(x)e-"«dx, (2.4) 

with the usual extension, by duality, to a mapping from S' to S' . 

In ( p.2[ ) the amplitude is written in this particular form to match the 

following definition: 



Definition 2.2. The energy- density of the solution of (2.1) is defined by 
n%x,t) :=\ip%x,t)\^. (2.5) 
We assume on the Weyl operator and on its symbol H: 
Assumption (Al) 

(Al)(i) 3 a eM.:H e S"^(M^) uniformly for e G (0,eo] 



(Al)(ii) 3 a unique self adjoint extension of iH^ {-,60) on 
By abuse of notation, we denote the unique s.a. extension by iH^ too. 



The hypothesis (Al)(i) means (see also po| ) that for all a,/3 E No, there 
exists Ca,/3 > , s.t. for all l,k G {l..m} and for all e G (0, eo] it holds 

\—-—Hix,0\ < Co^,pil+ I C \r-^ y{x,0 G K'^ X M'^. 

dx'^dCi 

In particular this implies that the Sobolev space H''{R'^) lies in the do- 
main of the operator . 

Furthermore this type of hypothesis are made in order to extend the rule 
of composition from differential operators to pseudo-differential opera- 
tors modulo lower order terms in e (for details see e.g. p3o[). 



Remark. We remark that the regularity assumptions on the symbol H 
are largely used for convenience and taken from |GMMP] in order to use 
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certain results, which were estabhshed there. They can be significantly 
weakened, for example: li H = |^P/2 + V{x) with V bounded below and 
V £ C^'^ (c f [ |LiPa[ ]), all results using Wigner transforms remain valid 
(see paMall ). 

On the initial data we impose: 
Assumption (A2) 

ni E L\R'^) ,ni>0 a.e. on M'^ and Sj G W/„'J(M'^). (2.6) 

Note that due to the low regularity assumed in (A2) a traditional WKB- 
expansion method would not be possible here ! This is one of the advan- 
tages of the Wigner formalism. 

Lemma 2.1. Let tl^j satisfy (A2) and assume (Al), then there exist a 
unique mild solution il)^{t) G C(Mt; L^(M'^)) of l\2. and its energy- 
density satisfies 

||n^(t)||^ = llnflli VtGM, (2.7) 
where \\ ■ ||i denotes the norm. 

Proof. Having in mind (Al)(ii) t he ass ertion is a simple consequence of 
Stone's famous theorem (see e.g. | ReSi[| ). □ 

Some particular examples for equation (^) are: 
Examples: 

(i) The Schrddinger equation 

edtij' -i—AiP^ + iV{x)7p^ = 0, (2.8) 

where e ~ ^. Here H{x, ^) = |^p/2 + V{x). 

(ii) The 1-d Airy equation (or linearized KdV equation) 

edti'' - y VLx = 0, (2.9) 
where H{x, ^) = Here e denotes the "physical" dispersion-parameter. 

(iii) The spinless Bethe-Salpeter equation (or " relativistic Schrodinger 
equation" ) 



£2 



edf^' -i\\l -—A + 1 + V{x) = 0. (2.10) 
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Again e ~ ^ and we have H{x,^) = y^|^P/2 + 1 + V{x). Note that in 
this example is a "true" pseudo-differential operator (i.e. the Weyl- 
symbol is not polynomial in ^). 

(iv) Another example of a "true" pseudo-differential equation assumes 
Hamiltonians of the form H = a(x)|^|, i.e. we have (in Weyl-quantized 
form) 

edti^^ + ie\D,M^Wiy))\^^,^ = 0. (2.11) 

In the constant coefficient case o(x) = 1, equations of this type can be 
traced back to the wave equation uf^ — Au^ = by noting that, the 
quantities 

^Pl{^,t) = dtu^±i\D\u'' 

satisfy 



We are now interested in the high-frequency limit e — > 0. For the sake of 
completeness we briefly review the traditional (naive) WKB-method in 
the next subsection. 



2.2 The WKB-method 

As stated in the introduction above we make the following ansatz 



ip^x, t) ~ A'ix, t) exp ( -S{x, t] 



X e 



(2.12) 



with > assuming (for the moment) that the phase and the amplitude 
are sufficiently smooth, and we expand the amplitude in powers of e: 

~ A + eAi + £^^2 + . . . 

We sketch the (formal) WKB-method for the case of a polynomial Weyl- 
symbol H = H(x, ^) with C°°-coefficients, i.e. 



H'^{x,eDMx)= J2^^'^D'y{ak{ 

\k\=o ^ 



x + y 



(2.13) 



y=x 



Substituting the representation ( 2.12 ) into (2J.) and collecting terms ap- 
propriately gives 



|fc|=o 



(2.14) 



Wigner vs. WKB 



8 



where Rk acts on as a differential operator of order \k\ < m (with 
coefficients depending on derivatives of S, H) . In particular 

Ro[A'] = H{x,V,S)A', 

^ ^ 90 dx,^2^^^dx,dxk d^.dCk 
" d^H{x,V,S) 



e 



-y 



where here and in the sequel we denote by the gradient w.r.t. the 
position variable, i.e. we consider y as a placeholder for the position 
variable x: S = S{y,t), H = H{y,(^) etc.. The last term in the expression 
Ri is obtained due to the fact that we use Weyl-quantized operators. 
Note that in the equations above all partial derivatives of the symbol H 
w.r.t. ^ a re evaluated at ^ = VxS. Plugging the above computations 
into ( p.l3| ), separating real and imaginary parts terms we obtain in the 
lowest orders 

— + Hix,V,S) = 0, 

dA dH{x,V^S) OA 1 A d'^S d'^H{x,V^S) 
dt dCj dxj 2 ^^dxjdxk dijdik 



1 

4e 



Note that the second equation is linear in A, it is called the transport 
equation for the amplitude. In terms of {n, S) = (A^jS) we obtain the 
following compact form, in the sequel called WKB-system associated to 



equations of type ( ^.l]) 



dtn + dw{nVi:H{x,V^S)) = 0, xGlR'^,tGM (2.15) 
dtS + H{x,V^S) = 0, (2.16) 
where the initial data are induced by the initial condition ( |2.2D 

n{x, 0) = 727(2;), S{x,0) = Si{x), xeR'^. (2.17) 

The firs t equ ation ( |2.15| ) is a conservation law for the energy-density, the 
second (|2.16| ) a Hamilton- Jacobi equation for the phase. 

To illustrate the method, we derive the (lowest order) WKB-systems for 



some particular examples of (2.1). 
Examples: 

(i) WKB-system for the Schrodinger equation ( |2.8| ): 

dtn + div(nV^5) = 0, x G M"', t G R (2.18) 

dtS+^-^^ + V{x) = 0. (2.19) 
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(ii) WKB-system for the 1-d Airy equation (2.9): 

dtn + d^{n{d^S)'^) = 0, x,t€R (2.20) 
a^S+lMl! = 0. (2.21) 

(iii) WKB-system for the (spinless) Bethe-Salpeter equation ( p. 10 ): 
f)^n + Hiv I _ ^^^"^ I =0, x£R'^,t£R (2.22) 



dtS+^^-^^ + l + V{x) = 0. (2.23) 
(iv) WKB-system for equations of type ( p. 11] ) : 



dtn + div ^ = 0, x G M'', t G M (2.24) 

dtS + a{x)\V^S\ = 0. (2.25) 

The equation ( p.25[ ) is the time-dependent eikonal equation of geometrical 
optics. Note that in one spatial dimension, i.e. d = 1, the above system 
de-couples and simphfies to: 

dtn + d^{nasgn{d^S)) = 0, x,teR (2.26) 
dtS + a{x)sgn{d^S)d,S = 0. (2.27) 

In this case the WKB-approximation is exact if dxS does not change 
sign, since then no terms ~ O(e^) appear in the expansion when using 
^fc = for A; > 1. 

<> 

In general, solving the system (|2.15 ), ( 2.16| ) allows an asymptotic descrip- 
tion of the solution of ( |2.1| ) in the pre-breaking regime, more precisely 
we have the following proposition, the proof of which is classic, see e.g. 
Q: 



Proposition 2.1. Letil)'^{t) solve (2.1), l\2.^ ), where is of the form 
{ 2.1c ). Assume ^Jrvf G 5(IR'*) and that Sj is bounded, together with all its 
derivatives up to sufficient order. Let < < tc^ be such that a smooth 
solution (n, S) of W^ ), §^ exists on M°' x {t^ , ) ^^"^ define 

^^fc,(x, t) := Vn{x,t)exp (^^S{x, t)^ . (2.28) 

Then, for every finite time-interval [ti,i2] C (icu^ca) '^^ have 

sup \\r{t)-r^,,m2<o{e), (2.29) 

ti<t<t2 

where 11 • lU denotes the norm. 
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Remark. The proposition can be extended to (systems of ) pseudo- 
differential operators, with somewhat weaker conditions on n, S and H 
(see e.g. pasl| , pas2| ). 



Generally one faces the following problem in this approach: 

The solution S of the Hamilton-Jacobi equation (l2.16| ) is obtained, at 
least locally, by the method of characteristics. This means that, for a fixed 
t, each point {y, t) G W in a (sufficiently small) neighborhood U C x Mt 
of the initial manifold x {t = 0}, is reached by a unique integral curve 
x{t,x), called ray, which can be found as described in the following. 
Let us define 

i{t, x) := VyS{x{t, x),t). (2.30) 



It is well known (see for example |lEv[| ), that the curves x{t,x),^(t,x) 
solve the IVP 

dx 

— = Vi,H{x,0, x(0,x)=x (2.31) 



-VyHix,i), iiO,x)=V,Siix). (2.32) 



dC 
dt 

Having solved this system we obtain the solution of ( 2.16| ) by integrating 

dS{x, t) 



dt 



V^H{x,i)-^-H{x,i), 5(x(0,x),0) =57(x). (2.33) 



For details see again [Ev|, [Fed|. As it is indicated above this theory is 
local, since it only gives the unique solution of ( |2.16| ) in a neighborhood 
lA of the initial manifold x {t = 0}. In other words, if we consider for 
every fixed t S M the map 

x{t, M'^ 

X ^ x{t,x) (2.34) 

it is, in general, not one-to-one for large times i G M, see e.g. the examples 
1.2, 1.3 in section 5. Further it is known that the phase is discontinuous 
at caustics, i.e. points of intersection of rays, which in general will hap- 
pen in finite times t = tc^^, t = tc^, with < < tea; called break-times. 
It is clear however that the formal WKB-expansion method can only he 
justified for tc^ < t < tc^ and thus a global, i.e. for all t € M, asymptotic 
description can not be obtained from this method. 



To overcome this deficiency, several generalizations of the method have 
been developed in the past decades. Most of them rely on the use of 
global Fourier- Integral operators (FIO) or the Maslov Canonical operator 
(MCO), which are now rather standard methods and covered in many 
books, however for the sake of the reader which may be not familiar with 
it, we shall now give a flavour of these methods. 
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2.3 The Fourier Integral Operator (FIO) 

Let us again start with an illustrative example: Consider the free Schrodinger 
equation in M'^, i.e. 

edtip' - i—AiP^ = 0, (2.35) 



^^(t = 0,x) = v/raK^e*'^^(^)/^ (2.36) 

Straightforward Fourier analysis shows, that its solution is explicitly 
given by the following oscillatory integral 

'>P%x,t) = {2TTe)-'^ [ / y^^e^^(^'^'«'*)/^dzde, (2.37) 
Jr Jr 

where 

S{x, z, t):={x-z)-C + ^lei' + Si{z). (2.38) 

Let us recall the theorem of stationary phase ([ |GuSt |, |lVIasl], |Mas2], 
ieS|): 



Proposition 2.2. Let a E C^(M'^), $ E C°°(M°') and assume that the 
set {z : 'V^{z) = 0, z £ supp(a)} consists of finitely many points zi, with 
i = 1 . . . N . If the Hessians D^{zi) are nonsingular, then for e <C 1 

(2^e)-^ [ aiz)e''^^'^^^dz = 
Jr 

, exp f -^(zi) + —mi I (a(zi) + 0(e)), 

where Zi = Zi{x,t) and mi := sgn(Z)2<I>(zj)) is the so called Maslov index 
of the the i-th ray. 

This implies that (locally) the main contribution to the solution of the 
Schrodinger equation stems from stationary points w.r.t. y and ^, i.e. 
points at which V^^^S = 0, which gives 

^ = V5/(z), x = z + tC, (2.39) 

i.e. we get a ray-map, c.f. ( |2.34 ), defined by the following relation 

X = x{t,z) = z + tVSi{z). (2.40) 

For small t the map z x{t, z) is singlevalued. In general however, there 
exist (maybe infinitely) many Zi = Zi{x,t), which obey the equation 
( 2.401 ). Note that the functions Sj{zi{x,t)) are local solutions of the 
Hamilton Jacobi equation 

O^S + ^-^^ = 0. (2.41) 
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Provided that the Hessians D'^S{zi) are nonsingular and assuming that 
there exist only finitely many points y^, the stationary phase theorem 
then gives us, the following multivalued WKB- approximation of if:^ 



N 



i=l 




nj{zi{x,t)) ( a { t \\ 



^^(x, t)c^y ^ W |, , : ' exp -Si{z,{x, t)) + -m, + 0(8 



tD^Si{z,ix,t))\ ^ ' ' 4 



Note that the mj change sign, each time a ray passes the caustic, which 
is defined here to be a point at which D^S is singular. 
For equations with variable coefficients, the above concepts have to be 
generalized, leading to the definition of a FIO acting on the initial datum 
V'f(x) by 

V''(x,t) ~ (27r£)-'^ / / A{x,z,^,t)e'^^'''''^''^/'il;j{v)dzd^. 



Since for t = the FIO has to reduce to the identity operator, we obtain 
the conditions 

S{x,z,^,t)\t=o = {x-z)-^, A{x,z,tt)\t=o = l- (2.42) 

At the core of the use of FIOs is the fact that the composition and 
commutators of two such FIOs is again "locally" a FIO (Hormander's 
Theorem [|Ho| ]) and that the rules to calculate their principial symbols 
are extensions of the rules of (composition and commutation of) pseudo 
differential operators. 

The leading contribution again stems from the stationary points, which 
correspond to the asymptotic solutions of the geometrical optics limit 
e — > 0. For points on the caustic the FIO is an integral which can be 
brought into a canonical form and evaluated in terms of special integral 
functions. This is the problem of classification of Lagrangian singulari- 
ties in contact geometry (see e.g. [ lAi| , [ AVG| ] ) . 



For the detailed implementation of these basic ideas we refer e.g. to | pu| ] , 
puSt| , |H§, 0, 0, 0, p^ ], p^ ]), however one should 
note that these approach generically assumes Sj,nj S C^(M'^), which is 
of course significantly stronger than our assumption (A2). 

After reviewing these classical methods of high frequency approxima- 
tion, we now turn to a new concept, namely the Wigner transformation 
technique. 



3 The Wigner function approach 

Definition 3.1. For given f,g(z S'{M.'^) and given e G (0, Eq] we define 
the Wigner-transform on the scale e by 

w^{f,g){x,0 = -^J f{x-e'-)-g{x + e'-)e'^^dz. (3.1) 
(2vr)'* J^d 2 2 
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For fixed e the phase space- valued function is a continuous, bilinear 
mapping: 

: S'{Ri) X S'{Ri) S'{Ri x Rf) 

In the following we shall use the notation: w^lf] = w^{f, /). 

The expression w'^[f] is real valued and for / G L^(R^) we conclude 

||u;^[/]||2 = (4vre2)-'^/2||p||2, (3.2) 

where the density "matrix" p is defined by p{x,y) := f{x)f{y). 

Remark. The Wigner-transform was originally introduced by E. Wigner 
in 1932 [Wi] in the context of semi-classical Quantum Mechanics. 
We remark that there are slightly different definitions of the Wigner 
transform in the literature depending essentially on the definition and 
normalization of the underlying Fourier transform. 

The inverse Wigner transform is given by 

p{x,y) ^ fix)f{y) = {27:)-'' j^w'if] {^,edj e^^^^-^^d^ 

(3.3) 

Note that this transformation allows to obtain the function / only up to 
a constant phase factor! 

Now a simple calculation shows that formally 

w'[f]{x,i,t)di = \f{x,t)\\ (3.4) 



which implies, that the energy density, as defined in ( |2.5| ), is the zeroth 
moment of w'^fV'*'] w.r.t. the velocity variable ^. This is sometimes called 
a microlocal decomposition of the energy provided by theWigner trans- 
form. A rigorous justification (Note that is in general not in L^.) is 
given by (see | LiPa | for details) 

n%x,t) = lim [ w'[i;']{x,tt)e-'^^^\'/^dC, (3.5) 

which converges in towards n"^ € L^(R'^). More generally, an impor- 
tant feature of the Wigner transform is that it facilitates a "classical" 
computation of expectation values (mean values) of physical observables 
{x,eD) in any state Tp"" , namely 

{i,',A^{x,eDW)^,= I [ A{x,Ow'm{x,OdxdC (3.6) 

(Here we assume that the symbol A{x, ^) is in 5(M'^'^).) In other words the 
real- valued function acts as an equivalent of the phase-space distribu- 
tion function, however in contrast to classical phase space distributions 
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the Wigner transform is in general not point-wise positivel Indeed it has 
been proved for example in [LiPa], |Hu] that w^[f] > if and only if, 
either f{x) = 0, or f{x) is a Gaussian. 

Remark. It is well known (see e.g. |LiPa|, |MaMa|, plu| ]) that averaging 
the Wigner function over phase space volumes large enough to fulfill the 
Heisenberg uncertainty principle yields a non-negative function. The so 
called Husimi functions are obtained by convoluting the Wigner function 
w^[f] in X and ^ with the Gaussian 

G'iz) := (7re)-2 exp f 

s.t. they become non-negative, i.e. Wj^{x,(,) := w^[f]-kxG^ ^^G"^ > a.e.. 

For the sake of illustration we present some particular examples of the 
"Wignerized" evolution equation (2.1): 

Examples: 

(i) The "Wignerized" Schrodinger equatio n (p.^) or Wigner equation (for 
details on this equation see e.g. | GaMa |, MaMa] ] ) 

dtw^ +^-V^w^ -e^[V]w^ = 0, x,CeM'^,tGM (3.7) 
where 9''[V] is the (non-local) pseudo-differential operator 

9^V]w%x,tt):=-^ [ [ [V{x + e'-)-V{x-el)] 

■w'{x,^,t) e'^^^-^'UzdC, (3.8) 



and w%t) := w^lip^t)], where ip^ solves (|j), (U). Note that in the 
free motion case, i.e. V{x) = 0, the Wigner equation becomes the free 
transport equation of classical statistical mechanics 

dtw^ + £, ■ V^w' = 0, x,(,£R'^;t£ M. 

Further note that in the case of potentials V which are quadratic in x the 
operator 6'^[V] simplifies to the classical force term for all positive 
e. The harmonic oscillator, V{x) = |xp/2, is a typical example 

dtw'' + ^ • Vxw'' - X ■ V^w" = 0, x,^eR'^,te M. 



(ii) The "Wignerized" 1-d Airy equation (2.9) 

dtw' + (^d^w' - —wlxx = 0, X, e, t e M. (3.9) 

This equation is obtained after a straightforward, but lengthy calcula- 
tion. 

We shall now perform the e ^ limit associated to the dispersive equa- 
tion (|0) with initial data (O), assuming (Al), (A2). 
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3.1 The WKB-limit using Wigner transforms 

One of the crucial properties (for a proof see [ GMMPf ) of the Wigner 
transform is 

Proposition 3.1. [ pMMP| ] If f and g lie in a bounded subset B of 
L^(M^), then w'^{f,g) is bounded uniformly in 5'(M^ x 
More precisely we have 

i^^^,w%f,g))ix,z) G CoiRi;L\R'')) 
uniformly for f,g £ B as e ^ 0. 



as e 



(3.10) 



Thus, by compactness, there exists a sub-sequence and a distribution 



w'^ E S' 



such that 



■ c/ 

w m 6 



(3.11) 



It has been shown (for example in | LiPa| ] , [ MaMa| ] ) that the limiting points 
of the Husimi function are also the limiting points of the Wigner func- 
tion. Since Wjj is non-negative, this implies that also is non-negative. 
More precisely it is a positive Bor el-measure, i.e. E Ai^, (the cone 
of bounded positive Borel measures) and therefore can be interpreted 
indeed as a c la ssical ph ase-space measure, called the Wigner measure of 
f'" (See |Gi, IGMMPI] for a proof of the positivity of u-^.). 

Remark. The Wigner measures are a particular version of defect 
measures, related to the H measures of L. Tartar [Ta] an d P. Gerard 



For more details see the expository article of N. Burq | Bu 



For simplicity we now restrict the scale e to a sub-sequence, such that 
the Wigner measure is unique, i.e. independent of the choice of the 
subsequence, and in the following we denote it by = w''. 
We further introduce the d-dimensional Poisson-bracket of two functions 
f,9, i-e. 



{f,9} := V,/ • - Vg/ • V,5. 



(3.12) 



In proposition 3.2 below we state the well known fact that the Wigner 
transform translates the action of an Weyl-operator asymptotically in 
zeroth order into a multiplication and in first order into a Poisson bracket. 
(A proof can be found in the appendix of | GMMP[ | ) 



Proposition 3.2. [pMMP|] Let P e for some a>0. Then iff,g lie 
in a bounded subset B o/L^(M'^), the expansion 



w'iP'^{x,eD)f,g) 



Pw'if,g) + -{P,nj^-{f,g)} + 0(e2) (3.13) 
2i 



holds in S'( 



uniformly for f,g E B. 
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Now let V'^ satisfy the IVP (2.1), ( |2.2| ). In view of ( |2.7D and proposition 
3.1 one obtains the uniform boundness of w'^[Tp'^]{t) in L°°(M; 5'(Mf x M|)) 

and thus the existence of w £ L°°(M; A^"'"(I 



'J^ X M|)) such that 



w^iip'] ^ w{dx,dC,t) in L°°(M;7W+(M^ 



weak-*. 



(3.14) 



Remark. It has been shown in |Gq| , |GMMP|| , |LiPa| that the hmiting 
process is actually locally uniform in t, i.e. w £ Cfe(K; M'^{ 

To derive an evolution equation for the limiting phase space distribution 
w we differentiate 

d^w' = w^ipt,^') +10^11;', ijf) = 2 Re '«;^(#,'0^). 



Now using equation (2A) and proposition 3.2, having in mind that 
and H{x, £,) are real-valued, we obtain the following essential equation 



dtw^ 



-{H,w'} + 0{e) 



(3.15) 



which is a linear PDE in w^. Passing to the limit e ^ yields, in the 
sense of distributions, the classical Liouville equation 



dtw + {H, w} = 0, X e M"', t G 



(3.16) 



We calculate the corresponding initial data wj which is the limit of the 
Wigner transform of the WKB-initial data ipf := ^/nje^^'^^ 



(27r)« 



ni{x - e^)ni{x + e|) 



exp 



J (^Si{x - e^) - Si{x + izi 



dz. 



(3.17) 



Lemma 3.1. Let w^[ipj] be given by iTf ), then 



wi := nj{x)6{£, — V xSi{x))dx. 



(3.18) 



Proof. The proof is a straightforward computation having in mind that 
the amplitude n/(x) and the phase Si{x) are by assumption e-independent. 

□ 

Remarks. 

• This special type of initial phase space distributions (|3.18D are called 
mono-kinetic initial data, i.e. for every x G there is exactly one 
(characteristic) speed vi{x) := VSi{x). Note that we have chosen 
the scale of oscillations in th e in itial data to be equal to the small 
parameter e in the equation ( |2.lD . 

We further remark that in case Si{x) G C°° the set {{x.,V xSi{x)y\ 
is a Lagrangian submanifold of phase-space M^'^, i.e. a manifold on 
which the symplectic form uj := dx Ad(^ vanishes (for details see e.g. 
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• The Wigner measure approach works for much more general initial 
data (cf pMMP[ , in this work we restrict to WKB data for the 
sake of comparison. 

Further we have: 



Lemma 3.2. ||GMMP|| Let w^[^lJ{t)] b e the W igner transform of the so- 
lution of (2.1) ^2. 2J ( calculated from Ii3.1(\ ), ^3.1$ )) and assume (Al), 



(A2), then the associated density n'^{x,t) given by ( \3.4 ) satisfies 



n^{x, t) n{x, t) :- 



w{x,dC,t) G a(IR'^;7W + (M^)) (3.19) 



where the convergence is locally uniform in t. 



Proof. A proof can be found for example in [GaMa], [GMMP], [ LiPa ], 
Note that the technical assumptions of "e-os cilla to ry" and "compact at 
infinity" initial data, which are introduced in |Ge|, |GMMP], are fulfilled 
if one imposes (A2). □ 

Remark. One can show that the e — > limit of other observables , 
or more general expressions which are quadratic in ip^ with a specific 
growth in ^, can be computed es sent ially in the same way, i.e. the limit 
is given b y the rig h t hand si d e of (|3.6D with u)^ replaced by w. For details 
see again IGaMaj , HGMMPj , |LiPa| . 



Thus by using the Wigner transform we obtained a (semi-) classical 
phase-space description, which we shall analyse in the next section. 



4 Analysis of the Wigner measure 

This is the main part of the work, in which we will try to establish 
the precise relation between WKB-asymptotic solutions of (2A) and the 
Wigner measure that has been obtained in the last section. We start 
with the following definition: 

Definition 4.1. The Hamiltonian flow Ft associated to the Liouville 
equation ( p. 10^ is given by 

Ft{x,0 = {x{t,x,0,i{t,^,0) (4.1) 
where (x^^) solve the initial value problem 

(IT ~ 

— = V^Hix,0, x(0,x,O=x (4.2) 
d£ . ~ 

^ = -VyH{x,o, e(o,x,e) = e- (4.3) 

Remark. The above ODE's are usually referred to as Hamilton's equa- 
tions and the curves (i,^) are often called bicharacteristics. 



Wigner vs. WKB 



18 



For the sake of simphcity we shah assume in the sequeh 
Assumption (A3) 

The Hamiltonian flow Ft is a continuously differentiahle globally defined 
map for every t G M. 

Remark. The global existence of Ft is a priori not guaranteed for gener al 



H G S'' , except for a < 1. (For more details see for example [|Ai1| , | ReSi[| .) 
We further remark that well known situations for which (A3) is valid, like 
the case of a-elliptic Hamiltonians H > Ci|^|" — C2 with Ci, 6*2,0 > 
are covered by our theory. 

An illustrative example where the Hamiltonian flow is not global in time 
is given in Appendix 1. 

A straightforward calculation shows that the Hamiltonian function H{x, ^) 
is constant along the flow Ft, i.e. 

H{x{t,x,^),^{t,x,0) = H{x,^), V t E M,x,e G K'^. (4.4) 

Also, by the classical Theorem of Liouville (see e.g. |Ar| ), we have that 
Ft is volume preserving; i.e. its Jacobian satisfies 

^^J dixit,xO,i(t,x,0) \^^ VtGM,x,eGM^. (4.5) 

With the above definitions the method of characteristics guarantees that 
the unique (weak) solution of the Liouville equation ( 3.16 ) subject to the 
initial condition w{t = 0) = wj £ A4"'"(Mf x M|) satisfies 

ipix,Ow{dx,dC,t)= [ [ ip{Ftix,0) wi{dx,dO (4.6) 



for all test functions ip G ^^(M^ x M|). In other words the solution w 
of (p.l6|) remains constant along all bicharacteristic curves (x,^). Since 



( 3.16| ) is a linear equation and (A3) is assumed we know that the above 



solution w exists for all t G M. In contrast to the c ase o f a nonlinear first 
order PDF we have global solutions for equation ( ^.16 ); i.e. no caustics 



appear in phase spacel This is sometimes referred to as "unfolding" of 
caustics. 



In the next subsection we connect the WKB-system ( 2.15[ ), ( p. 16 ) with 



a special class of solutions to the Liouville-equation, which we shall link 
to the traditional WKB method in the pre-caustic regime. 

4.1 On mono-kinetic phase space distributions 

Taking into account the monokinetic form of the initial data, we observe 
that, with the use of equation (|3.18| ), identity ( ^.6]) reads 

/ / ip{x,S,) w{dx,d^,t) = / nj{x)(p{x{t,x),£,{t,x)) dx, (4.7) 
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for all ip e Cb(R'^ x M^) where the curves x), ^(t, x)) solve the IVP 



( 2.31 ), ( |2.32| ). Note that, since the initial data are given by ^ = V^S'/(x) 



the curves x{t, x) are the rays associated to the Hamilton-Jacobi equation 
dtS + H{x,V^S) = 0, S{x,0) = Si{x). 

The above reflects the fact that the (x, ^){t, x, ^) bicharacteristics in phase 
space are projected down to rays x{t,x) in position space plus an addi- 
tional curve (,{t,x), which is the gradient of the phase S along the rays. 
More precisely, we have 

(x, i)it, x) = (x, i){t, X, V,Si{x)). (4.8) 
The function S is often called the generating function of the flow Ft . 

We shall use the following lemma for the proof of theorem 4.1. 

Lemma 4.1. Assume (A1)-(A3) and let w{t) he a (weak) solution of 
( 3. It ) with mono-kinetic initial data wj ( 3.18J . Then V t G M 



<i){H{x,i)) w{dx,di,t) = (j){H{x,vi{x))) ni{dx) (4.9) 

for every real valued continuous and hounded from below function (j) for 
which the right hand side of (^7^) is finite. 

Proof At first take e C(M) bounded. Setting ip{x,^) = cl){H{x,^)) in 



(4.7) combined with the fact that H{x,^) is conserved proves the claim. 



For unbounded cf) a density argument gives the result. □ 

Remark. If the Hamiltonian is bounded from below and gauged such 
that i?(x,^) > the above lemma states the energy conservation prop- 
erty of the Liouville equation, when setting 0(-) = id on M'^. 

Definition 4.2. We define the following set of functions: 

A := {A G C7(M;M+) : /" [ X{H{x, 0) wi{dx,dC) < oo}. (4.10) 
Now the following theorem holds: 

Theorem 4.1. Assume (A1)-(A3 ) and let -oo <tci<0< tc^ < oo. 
(i) The unique (weak) solution of (|3.16| ) 

dtw + {H,w} = in P'(R^ x M| x (tci,tcj) 
is given hy 

w{x, ^, t) = n(x, t)6{C - v{x, t)), (4.11) 
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if and only if the pair {n,v){x,t), where n G C;,((tci , ^ca); -^^(1^'^)) O'l^d 
V G L°°{{tc,,tc^);L\W^;dn{t)y) is a solution of 



dt{na{v)) +d\Y{n(j{v)V^H{x,v)) +nV^a{v)VyH{ 



(4.12) 



in V'iW^ X (tci , tea)); for every a G C-'^(IR'^) for which there exists a A G A 
such that: 



l + X{H) 



G L= 



kl|Vgg| 

1 + A(i7) 



G L= 



, |Vg^||V,g| 
^' 1 + \{H) 



G L°°(M'^). 
(4.13) 



(ii) M oreover if H is such that a = 1 and a = Vi, i = 1 . . .d satisfy 
(4-1^), then the monokinetic solution (4-11) is uniquely determined by 
the following system of fluid-type equations 



dtu + div(nV^//(x, v)) 
dt{nv) + div(V^-ff(x, v) nv) + nVyH{x, v) 

in P'(M'^ X (icii^cs)) with initial data 

n{x,0) = ni{x), u(x,0) = vi{x). 



0, 




(4.14) 
(4.15) 



(4.16) 



The eq uatio n ([4.12| ) is called a generalized formulation of the system 
( 4.14| ), ( |4.15 ). It reflects the fact, that different moments of the monoki- 
netic Wigner measure can be chosen, which is sometimes useful for nu- 
merical purposes. 



Proof. Choose a test function (p = 'y{x,t)a{(,) G 'D{M^ x 
then the weak formulation of ( 3.16| ) reads: 

/ / [{cT{Odt + a{OV^H{x,0-VMx,t) 

-VyH{x, • V^cr(07(x, t)] w{dx, dC, dt) = 
by inserting w given by ( 4.11[ ) we obtain: 

/ [{a{v{x,t))dt + (j{v{x,t))V^H{x,v{x,t))-VMx,t) 

—VyH{x,v{x,t)) ■ V^a{v{x,t))^{x,t)] n{dx,dt) = 

which is exactly the weak formulation of ( [4.12| ). Now let ai G P(M|) be a 
sequence of test functions satisfying the conditions stated in the theorem 
and converging a.e. to cr G C^( Mg) as / — > oo. Because the test functions 

(7/ satisfy the conditions above ( |4.13 ), the dominated convergence theo- 
rem (with the use of lemma 4.1. and definition 4.2. implies the assertion 
(ii). Now we can choose successively a = 1 and a = Vi, i = l,...,d in 
equation ( [1.12D to obtain the system ( |4.14D , ( [4.15| ) . 

□ 
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Remark. Note that in case H is given by 

Hix,C)=^-^ + V{x), x,^eR'', (4.17) 

where V is bounded from belo w we ca n ch oose A(-) = id and the above 
system of fluid- type equations ( 4.14| ), ( [4.15 ) simplifies to the well known 



zero-temperature Euler- equations of rarefied gas dynamics 

dtn + diY{nv) = 0, x G E"', t G (t^ , tcj (4.18) 
dt{nv) + div{nv (g) v) + nV = 0, (4.19) 
n{x,0)=ni, v{x,0)=vi, (4.20) 

which is, for tc, < t < tc^, the e ^ hmit of the QHD system (pT^), ( [Tsl ) 
that has been mentioned in the introduction. 

We now establish the precise connection between the WKB-system ( p.l5| ) , 
(pl^ ) and the ab ove theorem in the following corollary. 

Corollary 4.1. Let (tcntc^) be as above; 

(i) Further let S € C^(]R'^x(tci,tc2)) be a smooth solution of the Hamilton- 
Jacobi equation ^2.1(\ ) 

dtS + H{x,V^S) = 0, S{x,0) = Si{x). 

Define v{x,t) := 'VxS{x,t), and let n he the unique solution of 

then (n, v) satisfies 1 ^.1^ , ( [^. and is a generalized solution in the 

sense defined above. 

(ii) Let {n,v) S C^(M°' x (tcD^ca)) be a smooth solution of ^.14\ ), 



If the initial velocity is given by vi{x) = V xSi{x) G C^iW^), then there 
exists a phase function S{x,t), unique up to a constant, which is a solu- 
tion of l \2.l6i) on the same time interval. In particular the velocity v is a 
gradient field and the solution of ^S.ldJ can he written in the form 

w{x,^,t) = n{x,t)5iC - VccSix,t)) G Cbiitc,,t,,);M+{Ri x M^). 



Proof. Differentiating (2.16) w.r.t. Xi and using the chain rule (which is 
rigorous because we are in the regime of classical solutions) yields 



d 



d 



dt^S + j^Hix, V,5) + V^H{x, V^S) ■ ^V,5 = 



dy. 

and, with Vi := -^S we get 



(4.21) 



dtv + {V^H{x, v) ■ Vx)v + VyH{x, v) = 0. 

We multiply by 'V^a^v) to obtain 

a{v)t + {Vi^H{x, v) ■ Vx)a{v) + V^a{v)VyH{x, v) = 0. 

Again multiplying this equation by a regular solution n of the continuity 
equation (4.14) yields ( 4.12| ). Hence choosing appropriate test functions 
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we have proved that {n,v) is a generahzed solution on (tcijica)- We 
proove claim (ii), by first eliminating n from equation (4.12). Using 
(4.11) we obtain on (tci, ics) the equation ( [4.21 ) above subject to v{x, 0) 



\7xSj{x). One checks that the characteristic ODE system for this PDE 
IS given by (p3l| ), (12321) . Thus we can identify v(x{t, x),t) = ^{x, t) and 



conclude the existence of a unique (up to a constant) function S{x, t) s.t. 
V = VxS. The rest of the proof is identical to the one of (i) read from 
bottom up. □ 



Note that the system ^.14 ), l U-^^ defined n{t)-a.e. and thus remains 



valid as long as the density is single-valued I In particular this holds in 
the examples 1.2 , 2 in section 5 below. 

4.2 Multiple phases 

We will now show, under some more stringent assumptions, that away 
from the caustic we can locally extend w{t) beyond tci 2 ^ sum over 
mono-kinetic distributions, which will lead to a generalization of the 
asymptotic expression of the solution (t) of ( ^.l| ) . 



Since the Wigner measure remains constant along the flow Ft, we can 
write for e Cb{Ri x M|) 

{w{t), if) = {wi{F^t),v) = (n/(x(-t, X, OW.AO), ^) (4.22) 

where (•,•) denotes the duality bracket in the sense of measures and 
further we have used the following definition: 

Definition 4.3. The mapping 

fx,t{-) -.W^ is defined by 

fxM) ■■= i{-t,x,0 - VxSi{x{-t,x,0)- (4.23) 
Above (x,i) are parameters as the notation indicates. 
Definition 4.4. We denote the nullset of fx,t{') by 

/C(x,0:={eGM':/.,t(O = 0} (4.24) 
and the corresponding functional determinant by 

DfxM) ■■= det (^^) • (4.25) 

In general we cannot hope for fx,t{') to be a diffeomorphism in the whole 
propagation domain, at least we can get local results if we impose: 

Assumption (A4) 

Let ^ C X Mf be such an open set and N ^ N be such an integer that 
the nullset ( ^.2^ ) can be written as 



N 



lC{x,t) = \J{vi{x,t)] for all {x,t)eU. (4.26) 
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This means that we only allow situations with a finite number N of (gra- 
dients of) phases at each point of the propagation domain. Of course w 
is well defined even if (A4) does not hold, however in general no upper 
bound N can be found, = oo is possible! 

Remark. Indeed very little is known about this problem in general 
situations, so far the only well studied example is the one dimensional 
free motion case, i.e. H = -ff(^), ^ G M. Under the additional assumption 
that H{^) is strictly convex, i.e. there exists c > such that H" {^) > c, 



Tadmor and Tassa have shown in [TaTa], that if H'(Sj{x)) has a finite 



number of inflection points, then the number number of "original" shocks, 
i.e. shocks that do not result from the interaction of other shocks, in the 
entropy solution to the corresponding conservation law is finite, which 
implies (A4). An example of an initial condition VSj which evolves, in 
the particular case of H = ^^/2, into a n a.e . C°° function with countably 



many original shocks can be found in | Sch | 



Theorem 4.2. Assume (Al)-(A^) and denote by Vi{x,t), i = 1...N, 
the elements of the nullset IC{x,t). 

If the point {x,t) £ U is such, that Dfx^t{vi) ^ for all i = 1...N, 



then the measure- valued solution of the Liouville equation (3.16) can be 
written as 

t t +\ x^ni{x{-t,x,VxSi{x,t))) 
H^,^,t) = Y^ |I^/..,(V.5.(x,t))| ^(^-^-^-(-'^)) (^-27) 

where Si £ Cf^jU) such that VxSi = Vi £ IC{x, t) for all i £ {1 . . . N}. 

Proof Choose tests functions ip{x,C,t) = j{x,t)a{S,) £ Cb(Mj; C6(M^'^)). 
As above we write 

{w{t), ^) = {ni{x{-t, X, OWxAO), ^) ■ (4.28) 

By the coarea formula (see [ pv| , p^]) we have 

{5{fxM))^<^m = {m,^.A0) 

with 

■■= [ ^ ^(Ottt^^, (4.29) 

which is well defined as long as Df^AO / 0) even if fx,t{0 is not an iso- 
morphism for all points £ By the definition of the delta distribution 
the above is equal to 



K.{x,t) 



\DfxAoy 



since {f~l{0)} = £ M.'^ : fx,t{0 = 0} =: }C{x,t). By assumption 
(A4) we have that /C(x, t) is a finite union of points Vi{x, t), thus we can 
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evaluate the integral in terms of 



/ 



N 



sr^ a{vi{x,t)) 



{S{C-v^{x,t)),a{0) 



which in the sense of measures can be written in the form 

N 

^ {0' 

The local differentiability of Vi, i = 1. . .N, is a direct consequence of 
the implicit function theorem. The existence of phase functions Si{x,t) 
such that Vi = VxSi can be concluded essentially in the same way as in 
Corollary 4.1, using 

fxAO = C - V,.S/(x) (4.31) 
and thus 

fx,oM^^ 0)) = Mx, 0) - V^Siix) = 0, (4.32) 

since by construction Vi{x,t) G }C{x,t). Thus we obtain ([4.27 ) by insert- 
ing (0^ ) into (HH). □ 

In view of the above theorem we can now define define the caustic set 
C C X Rt by 

Definition 4.5. 

C := {{x, t) : Dfx^tiO = for at least one i G lC{x, t)] (4.33) 
This definition becomes more clear by the following lemma. 

Lemma 4.2. Denote by x^^{t,x) C R'^ the pre-image of the point (x,t) 
under the flow x{t, •) (i.e. the set of all points z with x{t,z) = x) and 
by x{t, •) the (locally defined) inverse map of x{t, •). Let V C x M.t be 
such that there exists only one v = V^S G IC{x,t), then we have 

\Df.,t{VxS)\ = J{x{t, x),t), V (rc, t)eV (4.34) 

where 

J(z,t):=|det(^^)|. (4.35) 

Thus the Jacobian J of the rays x{t,x) becomes zero at the caustics, 
c.f. theorem 4.3. This property is usually used to define caustics in the 
WKB-framework (see e.g. [ pKM|[ , [ pfedf ). 



Proof. The claim follows directly from identity ( [4. 61 ) , having in mind that 
if there exist only one v G IC{x,t) then x{t, •) is locally one-to-one and 
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thus we can apply the transformation law of integrals (see pq]) on the 
right hand side of (4.6), to obtain 

n{x, t) = ni{x{t,x))\ det ' ' ^ 



dx 

Since x{t, •) in this case is, by assumption, invertible the basic calculus 
for determinants implies 

nr(x(t,x)) 

n{x,t) = V \ \ - (4-36) 

^ ' ^ J{x{t,x),t) ^ ' 

Finally note that f (t, x) = x{—t, x, v{x, t)) with v = VxS{t, x) as long as 
x{t, •) is one-to-one. Thus the claim is proven. □ 

Note that ( [4.36| ) is the solution of the conservation law ( p. 15 ) before the 
caustic onset, i.e. for tci < * < ^ca- 

We have obtained a multivalued description of the solution of the WKB- 
system for all t G R, locally away from caustics C, by representing the 
limiting phase space density w{t) as a sum over mono-kinetic terms, 
each of which can be associated to a sing le br anch of the multivalued so- 
lution of the Hamilton- Jacobi equation ( |2.16| ). Note however that, since 
VxSi = V xSi with Si{x,t) = Si{x,t) + Cj, Cj G M, it is clear that for 
each Vi £ IC we obtain the corresponding phase Si only up to a constant. 

Remark. Physically the multivaluedness can be interpreted as interfer- 
ence of the wave with itself, or in terms of classical mechanics, that faster 
particles overtake slower ones for times t < tc^^, t > tc^- 

To close the argument, the following corollary shows, that the Wigner 
measure given by ( |4.27 ) is indeed the limiting measure of a superposition 
of (WKB) waves. 

Corollary 4.2. Assume (A1)-(A4)- and letUnC = 0, with C defined in 
^.33[) . Define for all {x, t) £ U a generalized WKB -asymptotic solution 
by 

^tnkbix, t) ■■= J2 VMx,t)exp i^^Siix, t)j , (x, t)eU (4.37) 

where VxSi £ K., M i = 1 . . . N and each branch of the energy density 
reads 

I ,N nj {x{-t,x,VxSi{x,t))) 
"^^"'^^•= \DfxAy.S.ix,t))\ ■ ^'-''^ 



For a localization let (p G Cq°(M^. x M^) with supp (j) ^U. Then the u nique 
(semi-) classical Wigner measure '^f '^'^{4''4^wkb\ '^^ given by ( 4-27 ) mul- 
tiplied by (fP' . 
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Proof. Having in mind lemma 3.1, the proof follows directly from the fact 
(see [ p^ ) that if 5^ have Wigner measures Wf, Wg, which are mutually 
singular, i.e. there exist two disjoint Borel-sets A, B with Wf{A^) = 0, 
Wg{B'') = 0, then 

^e, ye ^ + 5'(M^ X M|). (4.39) 

Clearly this is true in our case, since each term in the expression ( [4.27| ) 
is a measure which is supported on {(x, Va;5j(x, t))}. By construction, 
we have VajS'j 7^ xSj for i ^ which implies the assertion. 

□ 

Observe that V'^fct does not explicitly contain the so called "Maslov-phase 
shifts" (c.f. section 2.3) due to the fact that each branch of the multival- 
ued phase Si, obtained by our kinetic approach, is uniquely determined 
only up to a constant, which is not explicitly specified by the approach. 
In each region U the Maslov indices are such constants. More precisely, 
we connect the notation used in our Wigner measure approach with the 
one used in section 2.3 (FIO's) by stressing that 

Zi{x,t) = x{-t,x,VxSi{x,t)), 

evr 

Si{x,t) = Si{zi{x,t)) + -^rui. 



On the other hand it is clear, see section 2.3 and e.g. |GuSt] for more 
details, that once the the Lagrangian manifold generated by V xSi{t,x) 
is given, the Maslov phase shift can be calculated from it from purely ge- 
ometrial considerations. Thus although it does not appear explicitly, the 
Wigner measure (in the representation ( 4.27| )) contains all information 



needed to obtain the phase shift, which then allows the construction of 
the correct multivalued WKB-expansion after the caustic. 

Remark. Although the Wigner measure gives the correct multival- 
ued description after the first caustic its computation in general is quite 
labourintensive, in particular since it involves 2d + 1 variables. From a 
numerical point of view one would like to work in physical space x . 
This requires the approximation of u) by a system of 2N equations for the 
2N unknowns {ni,Vi), 1 < i < N . This system for the 2N moments of w 
in general is not closed. However, in geometrical optics, the multivalued 
form of the Wigner measure (for a fixed A^) gives a closing condition 
which allows (in principle) the correct description of multivalued situa- 
tions until the next caustic forms. For details on this proble m (in d = 1) 
see | JiLi | and for some alternative approaches we refer to [Bel], | Be2| ], 



pa. 



4.3 Concentration effects 

We now describe the behavior of the density at focal points, which typi- 
cally arise as the onset of caustics. We will be able to distinguish between 
two specific cases of energy concentrations. 
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Theorem 4.3. Fix t ^ 0,y £ and let V ^ x Rt be a region 
(closed set), with {y,t) G V, such that there exists only one v G IC{x,t) 
for all {x,t) E V. Further assume that there exists r > such that 
{{x,t) : \y — x\ < r} C V. Then we have 

Hi {x{-t,x,v)) 

n{x, t) = ^^^^ ^^^^^^ ^^^^ X{x^y} + - y), \y-x\<r (4.40) 
where x denotes the characteristic function and 

Ijl:= I ni{z)dz. (4.41) 

Proof. Applying identity ( [4.6D to the monokinetic form of the initial data 
with ip{x,S,) = 7(x) we obtain 



j{x)n{dx,t) = I nj{x)'j{x{t,x))dx 

Thus we can write 

j{x)n{dx,t) = I nj{x)'y{x{t,x))x{x^£-i(y,t)}dx 



+ / ni{x)-f{x{t,x))x{xex-^{y,t)}dx. (4.42) 

JR'' 

Now let 7(x(t, •)) be supported in {x R'^ : \y — x\ < r}. In the last 
term of the right hand side we obtain 



7(y) / ni{z)dz = {fi6{x - y),j{x)) 

where (•,•) denotes the duality bracket. In the first term on the right 
hand side the transformation law of integ rals can be applied thus using 



lemma 4.2 we obtain the expression ( 4.40| ) which proves the claim. 

□ 

If there is a nonzero amount of initial mass /i carried by rays int o the 
point {y,t) the energy density n "concentrates" as the last term of ( [4.4[l| ) 
shows. This is also an explanation of the word "caustic", since its Greek 
origin means "which burns". We shall refer to these caustic points as 
"hot" 



We deduce from the theorem above the following easy consequences: 
Corollary 4.3. If wi{{x~^{t,y)} x M^) = 0, then locally around y 

, , ni(x(t,x)) / N /. . N 

^(^'*) = 7^|^X{.^.}' (^,i)eV. (4.43) 

In particular the density remains in LL( R'^) in this case. Clearly ifx{t, •) 
is a diffeomorphism in V the expression ^.45 ) is valid for all {x,t) G V. 
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Proof. The claim follows from nj{A) = iff wi{A x M^) = 0, where A is 
an arbitrary Borel set. □ 

It is important to note, that although the amplitude blows up at the 
caustic C, the density n, still makes sense as a measure. For this reason 



the amplitude is sometimes called half- density [Ho| 



In other words: Although some rays x{t,-) may cross, the density, as 
the corollary shows, may still be in Lj^^lW^). This will be referred to as 
a "cool" caustic point, in contrast to a hot caustic where we obtain a 
concentration of the density. A particular case of such a cool caustic is 



given in example 1.3. in section 5 below. (See also [IVIPP | for a numerical 



study.) Another cool caustic is given by example 1 in [paMa |. 



Remark. Clearly the above theorem and corollary can be extended 
to regions in which finitely many (gradients of) phases appear. This is 
in particularly true for points on the one dimensional cusp-caustic (see 
example 1.3), where within the caustic region we have 3 zeros Vi{x,t) 
and outside there is only one. Thus we obtain cool focus points on each 
branch of the caustic, whereas the caustic-onset point, or focus point, is 
hot. 



5 Case studies 

We now illustrate the above analysis with examples. 
1. Free motion 

One should note that the following calculations are generalizations of the 
formulas given in the example of the free Schrodinger equation in section 
2.3. 

Although the generalized free motion case, i.e. 

H = HiC), (5.1) 

is the most simple one, it nevertheless features a remarkable variety of 
interesting phenomena. The associated Hamiltonian flow is given by 

Ft{x,0 = {x + tV^H{^),0. (5.2) 

The velocity remains constant ^(t, x,^) = ^, Vt G M, x,^ G and 
thus the rays are straight lines, with slope 'V^H{\/xSi{x)), i.e. we have 
x{t,x) = X + tV {V xSi{x)) . In this particular case definition 4.2 reads 

f.AO ■■= e - V,S/(x - tV^HiO) (5.3) 
and its zeros Vi{x,t) £ IC{x,t) satisfy the implicit relation 

Vi{x, t) = VxSiix - tV^H{vi{x, t))). (5.4) 
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This is the well known (multivalued) implicit solution formula of the 
conservation law 

dtv + V^H{v) = 0, xeR'^,teR (5.5) 
vi{x,0) = V,Si{x), (5.6) 

which holds as long as the determinant of the Jacobian is nonzero 

Df.,M) = det(/ + tD^Siix - tV^H{vi))D^H{v,)) / (5.7) 

(Here D^f denotes the Hessian of /.) The multivalued phase Si{x,t) 
is obtained, locally in each region U in which is constant, by using 
standard Hamilton-Jacobi theory, i.e. by integrating (2.33) 

^^l^ = H'{i).i-H{i), S.{x,0) = Sjix) 

after inserting ^(t,a;) = Vi{x{t,x),t) and using Vi{x,0) = VxSj{x,0) to 
determine the initial condition. 

Finally, it follows from ( [4.27| ), that for all {x,t) G U, {x,t) C, the 
density is given by 

n(xt) = Y' ^dx-tViH{v,)) 

^ I det(l + tD^Siix - tV^H{vi))D^H{vi))\ ' ^ ^ 

In the examples 1.1.-1.3. below we restrict ourselves, for simplicity to 
the case of one spatial dimension, i.e. d = 1. We further assume that 
i?(^) is equal to the classical kinetic energy 

^(0 = f, eeiR, (5.9) 

which corresponds to case of the free Schrodinger equation. Then we 
analyse for t > the behavior of the WKB-system subject to different 
types of initial phases Sj. 

Example 1.1. 'Wo caustic" 

Sj{x):=^, xGM (5.10) 

Here the rays x{t,x) never cross, instead they spread out, forming a so 
called rarefaction wave{see fig. 1). The only element v{x,t) G IC{x,t) is 
given by 

The solution of the Hamilton-Jacobi equation, which is single valued and 
smooth for all times t > 0, reads 
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and the limiting density ( |5.8D simplifies to 

nix,t) = -^nj(-^]. (5.13) 



|t + 1| \t+l. 
Example 1.2. "Focusing at a point" 
By simply changing the sign of the initial phase we obtain 

2 

Si{x):=-^, xGR, (5.14) 

which leads to the single focus case. All rays intersect at the hot focus 
point {x,tc) = (0,1) and spread afterwards (see fig. 2), i.e. there is a.e. 
only one phase 

2 

5(:E,t) = ^^^, t^l. (5.15) 

The above theory does not tell us anything about the precise description 
of the phase S{x,t) at t = 1. With the use of the theorem 4.2 and 
equation (^) the density is given by 

n{x,t) = ^l^n/(^^) t/1, (5.16) 
n{x,l) = / nj{x)dx 6{x). (5-17) 



In this example there exists a generalized solution (n, (7{v)) of ( 4.12 ) n(t)- 
a.e.. 

Example 1.3. "Caustic" 
We choose 

Sj{x) := - ln(cosh(2;)), xGR (5.18) 

such that the initial data for the ODE-system of rays is " compressing" : 

^^(x) := -tanh(x). (5.19) 

The equation which characterizes the kernel IC{x,t) cannot be solved 
explicitly, however a precise numerical study is given in [ [MPP ]. We want 



to stress again that in this case no "hot" focus appears, i.e. the limiting 
density n G L}^^{R) for ah t e M. 

An explicitly solvable example (see also [ GaMa[ |) which has a similar 



qualitative behavior (except that the focus is hot) is given by 

Si{x) ■.= xx{x<o} + ix- —)X{o<x<i}, xGR (5.20) 
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where x denotes the characteristic function. Note that the initial condi- 
tion is only Lipschitz continuous, i.e. differentiable almost everywhere. 
Up to the time t = 1 the rays do not intersect, at i = 1 a focus point 
occurs at {x,tc) = (1,1) from which two caustics emanate, forming a 
cusp (see fig. 3). By applying theorem 4.2, in a neighborhood of the 
focus, we obtain for the density 

n{x, 1) = ni{x)e{x - 1) + ni{x - 1)G(1 - x) + ^i6{x - 1). (5.21) 

Here the amount of mass that gets concentrated at the point (1, 1) is 
given by // := ni{x)dx and Q{x) denotes the Heaviside function. 
For t > 1 the solution of the Hamilton-Jacobi equation is triple-valued 
(as it is generic for the singularities in one dimension) within the region 
1 < X <t, since fx,t{-) has three zeros there 

vi{x, t) = 1, V2{x, t) = V3{x, t) = (5.22) 

and thus we get for the density within the caustics 

1 / X — t\ 

n{x, t)x(i<x<t) = ni{x -t)+ _ ^^ nj [j—tj ^^'^^^ 

The corresponding phases are obtained by a simple integration using 
( 2.33|) , for example the phase function, that gets transported into the 
caustic region from the left, is given by 

Siix,t)=x-^. (5.24) 

We now turn to an explicitly solvable case with x-dependent Hamiltonian 
H{x,0- 

2. Harmonic Oscillator 

Consider in the example of the harmonic oscillator (see also | GaMa| ) 

H{x,0=^-^+^-f- x,eGE^ (5.25) 

In this case the Hamiltonian flow Ff is given by 

x{t,x,(,) = xcost + (^s'mt 
^(t,j;,^) = —xsint + S^cost. 
If we choose in particular 

Si{x):=kx, k>0,x£R (5.26) 

we obtain a constant initial velocity vj{x) = k > 0, such that all rays 
intersect at hot focal points (see fig. 4) given by 

{x,tc) = {{-ir^^k,{2m + l)TT/2) meZ. (5.27) 
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We obtain a.e. only one zero ^ = v{x,t) S fC{x,t) of fx,t{C)^ given by 

k — X sin t , , , vr 

v{x,t) = t/(2m + l)- 

COS ^ 

and thus the solution of the Hamilton Jacobi equation is unique n(t)-a.e. 

S(x,t) = --ix"^ + k^)tant + t/(2m + l)-. (5.28) 
2 cos t 2 

The density is given by 

1 /x — A; sin A , , , vr , , 

n{x,t) = —^ni( — , t/(2m + l)-, (5.29) 

|cosip \ cost / 2 

n(x,t) = / ni{z)dz 6{x - {-l)"'k), t = (2m + l)J. (5.30) 



In this case again there exists a generalized solution {n,a{v)) of ( 4.12| ) 
n(t)-a.e.. 



6 Conclusion 

In this paper we compare two important approaches within the theory of 
geometrical optics for linear dispersive operators (acting on namely 
the Wigner transform and time-dependent WKB-asymptotics. Whereas 
the latter faces the problem of caustics the limiting Wigner measure is 
insensitive for such obstacles, due two the fact that it lives in a 2d- 
dimensional phase space in which the appearing singularities become un- 
folded. This feature makes the Wigner transform the method of choice 
if one is particularly interested in the high frequency behavior of energy 
densities, which are obtained as moments (w.r.t. the velocity variable ^) 
of the Wigner measure. 

Maybe the most important conclusion from this paper is the fact that for 
WKB-initial data the limiting Wigner measure w can be (locally away 
from caustics) decomposed into a sum of monokinetic terms ( [4.27 ). Each 



term carries all the information needed, namely (n^, 5j), i = 1 . . . A^, to 
obtain a (local) approximation V'^^^ of the solution ijj'^ to the dispersive 
equation at the order 0(e), despite the fact that the measure w does 
not explicitly contain the phase shifts obtained by the stationary phase 
method used in the FIO approaches. However since it is well known that, 
given iSj, these phase shifts can be obtained by a purely geometrical com- 
putation, we conclude that these important feature of geometrical optics 
is actually hidden in w. 

In other words the Wigner measure which is a description for e = does 
indeed allow to obtain a approximation of ip'^ for finite e and can thus be 
seen as a necessary (but of course not sufficient) condition for the validity 
of the multivalued WKB approximation. Higher order correction can be 
also obtained from the Wigner transform (using the expansion ( ^.13 )), 
although we neglected them in this work for the sake of simplicity. 
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We further stress the fact that although the amphtude of the WKB ap- 
proach becomes infinite at the caustic, the energy density n ("square of 
the amphtude" ) remains well defined in the sense of measures, as theorem 
4.3 shows. What can not be obtained from the Wigner measure is infor- 
mation on the phase at caustic points, which is clearly in agreement with 
the fact that the multivalued WKB approximation breaks down there 
too and the fact that the only correct description at caustics is given by 
FIO's. 

It remains to say that the Wigner transform can be used in situation 
with much less regularity of the initial data than the WKB approach or 
its generalizations, since all our computations remain valid for nj G 
and Sj G with Lipshitz continuous derivative. 
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7 Appendix 

As an appendix to the main results of this paper we first present a (rather 
exotic) example where key assumptions of the presented theory are vi- 
olated. Then we briefly comment on the physical interpretation of the 
fluid type equations, which arise in our work. 



7.1 Appendix 1 : A counterexample to global Hamilto- 
nian flow 

As a "caveat" that the assumption of essential self adjointness (Al) (ii) 
and global Hamiltonian flow (A3) are not trivial we consider a variable 
coefficient Airy-type equation ( p. 91 ) in one space dimension 

H{x,0 = -xf, x,^eR. (7.1) 

A lengthy calculation shows that the Weyl-quantized Hamiltonian oper- 
ator corresponding to this symbol does not have a self adjoint extension. 
Moreover the Hamiltonian flow Ft associated to (7.1) is only locally de- 
fined as an explicit calculation shows: 

i(t.i,0 = i:(l-2ftf''\ (7.2) 



We remark that to our knowledge it is an open question if essential self 
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adjointness of {.,eD) is a sufficient condition for global Hamiltonian 
flow defined by H{x,S,). 
If we again choose 

Siix) ■.= kx, k>0 , xeR, (7.4) 

all rays x{t, x) focus at the point (x, tc) = (0, 1/2A;^), i.e. we again obtain 
a single focus case as in example 1.3. above. The corresponding phase- 
function is given by 

S(x,t) = t<^. (7.5) 

It is not possible to extend the solution beyond the break time since 
Ftixj^, = k) is only defined for t < tc- 

In this example neither the WKB nor the Wigner approach give an 



asymptotic description of solutions of (2J), (|2.2| ) after breaktimes. 

We further remark that changing the sign in the Hamiltonian function 

gives rise to an expansive flow of rarefaction-wave type. 

7.2 Appendix 2: On the arising fluid-type equations 

As indicated above the fluid- type system ( 4.14 ), (|4.15| ) admits a physical 
interpretation if the Hamiltonian function is given hy H = |^P/2 -|- V{x) 
(see the remark below theorem 4.1). In general this is not the case. 
However we state in the lemma below that one can find indeed an equiv- 
alent system which is of the same form as the Euler-equation of gas 
dynamics, and which can be defined for general Hamiltonian functions. 

We define the generalized velocity u and the (time dependent) modified 
force term / by 

u{x, t) := V^H{x, v{x, t)), (7.6) 

where v{x, t) is a d-dimensional vector field. With this definitions we 
obtain the following result: 

Lemma 7.1. Let {n,v) be a smooth solution of ( J^.l^ ), o,f^d let 



u, f be defined as above on the same time interval (tc-^jtc^), then {n,u) 
satisfies 

dtn + div{nu) = 0, n{x,0) = nj{x) (7.8) 

dt{nu) + div{nu u) + nf = 0, u(x,0)=u/(x) (7.9) 

with ui{x) := V ^H{x^vi{x)). 

Proof. By definition (|7.6D, it is clear that the conservation law (7.5) is 



equi valent t o (4.141 ). Further note, that after eliminating n{x,t) from 



(|7.9D using (7.8) one obtains for u{x,t) the Burgers equation with source 
term 

dtu + iu-Voc)u + f = 0. (7.10) 
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Thus it remains to show that if u{x,t) solves (7.10), then v{x,t) is a 
solution of 



dtv + {V^H{x, v) ■ V.^)v + V^H{x, v) = 0, 



(7.11) 



which is again obtained from ( |4.15p by elimination of n{x, t). Calculating 
inner derivatives we obtain for the i-th component of ( |7.1C| ) 



k=l 



1=1 



1=1 



In order to have a classical solution v{x,t) of ( [7.11| ) this implies that 
each component , i = 1 . . . d, of the modified force term must satisfy 
the relation 

d d 
k=l 1=1 

which is exactly definition (|7.7|) above and the claim is proved. □ 



Clearly the fluid-type system formulated in (n, v) and the Euler system 
in {n,u) are equivalent (for smooth solutions) if the relation (7.6) can be 
uniquely solved for v in terms of li. In case 

H = u;{0 + V{x), (7.12) 

where a; is a general dispersion relation, we obtain u = V^uj{v) and 
/ = D'^uj{v)V^V (where again D'^uj denotes the Hessian matrix of uj) 
and thus 

dtn + div{nu) = 0, n(x,0) = rij (7.13) 
dt{nu) + div(nn u) + nD'^uj{v)VxV = 0, u{x, 0) = ui. (7.14) 

Further note that by definition u in general is not the gradient of a phase 
function S. 
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